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In  1939  Weibull  applied  weakest  link  theory  to  the  interpretation  of  the 
variability  of  fracture  stress  of  nominally  identical  specimens.  He 
attributed  this  variation  to  the  presence  of  sources  of  weakness,  or  flaws 
that  were  given  a strength,  but  no  other  physical  properties.  Accordingly 
bis  theory  rests  primarily  on  the  statistics  of  extreme  values. 


ir  is  now  generally  believed  that  the  most  important  sources  of  weakness 
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are  microcracks.  Thus  an  ideal  statistical  theory  of  fracture  must  take 
into  account  of  fracture  mechanics  and  material  microstructure.  Present 
paper  outlines  recent  progress  along  these  lines,  and  shows  how  the  re- 
sulting more  fundamental  theories  reduce  to  Weibull  theory  when  certain 
simplifying  assumptions  are  made. 
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ABSTRACT 

In  1939  Weibull  applied  weakest  link  theory  to  the  interpretation  of 
the  variability  of  fracture  stress  of  nominally  identical  specimens.  He 
attributed  this  variation  to  the  presence  of  sources  of  weakness,  or 
flaws,  that  were  given  a strength,  but  no  other  physical  properties. 
Accordingly,  his  theory  rests  primarily  on  the  statistics  of  extreme  values. 

It  is  now  generally  believed  that  the  most  important  sources  of 
weakness  are  microcracks.  Thus  an  ideal  statistical  theory  of  fracture 
must  take  into  account  of  fracture  mechanics  and  material  microstructure. 
Present  paper  outlines  recent  progress  along  these  lines,  and  shows  how  the 
resulting  more  fundamental  theories  reduce  to  Weibull  theory  when  certain 
simplifing  assumptions  are  made. 


INI  RODUCTiON 

In  1939,  * ^ Weibull  intioduced  weakest  link  theory  into  mechanics  with 
the  primary  purpose  of  accounting  for  experimentally  observed  variations 
in  the  fracture  stress  of  nominally  identical  specimens.  He  showed  how 
test  data  on  laboratory  specimens  uniformly  stressed  in  simple  tension 
can  be  used  to  find  the  failure  statistics  of  specimens  of  different  size 
and  nonuniform  loading.  By  assuming  that  only  tension  contributes  to 
fracture,  he  extended  the  theory  to  apply  also  to  bending,  and  in  fact  to 
any  uniaxial  combination  of  bending  and  tension.  Because  of  the  essen- 
tially rational  nature  of  the  theory  and  its  generally  good  agreement  with 
experiment,  it  has  become  wTdely  known  and  used.  For  convenience  this 
body  of  theory  will  be  referred  to  herein  as  Theory  A. 

Weibull  also  gave  a procedure  for  finding  the  probability  of  failure  for 
finding  the  probability  of  failure  for  polyaxial  stress  states.  This  procedure, 
which  we  shall  call  Theory  B,  requires  a treatment  of  uniaxial  tension  somewhat 
different  from  that  employed  in  Theory  A.  For  instance,  the  number  of  weak 
spots  per  unit  volume  found  using  Theory  B differs  from  that  found  using 
Theory  A. 
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Theory  B has  been  questioned  by  a number  of  investigators.  Barnett,  et  al. 
considered  Weibull1 s rule  for  polyaxial  tension  to  be  plausible  but  not 

necessarily  correct,  and  listed  five  similar  procedures  they  considered 
equally  plausible.  Some  investigators  who  accept  Theory  A choose  not 
to  use  Theory  B,  and  employ  some  other  approach  to  polyaxial  stress 
states,  such  as  the  assumption  of  independence  of  principal  stresses. 

It  may  be  that  Weibull  came  to  doubt  Theory  B himself  since  he  wrote 
in  1966,^  "Another  problem  of  a more  theoretical  nature  will  be  to  deduce 
the-  effect  of  bi-  and  triaxiality  on  the  distribution  functions  of  one- 
dimensional stresses.  If  the  principal  stresses  are  acting  independently 
of  each  other,  and  it  seems  that  such  materials  may  exist,  then  the 
problem  may  be  soluble  along  lines  previously  sketched.  In  other  cases 
the  solution  is  very  intricate  and  will  certainly  require  close  examination 
of  the  physical  behavior  of  the  material  in  question.  " 


Weibull  assumed  that  the  fundamental  cause  of  the  variations  in 
fracture  stress  was  the  presence  of  weak  spots  or  flaws  in  the  material. 


-2- 


He  suggested  that  they  might  be  due  to  minute  foreign  bodies,  impurities, 
or  cracks.  The  flaws  were  given  strengths,  but  no  other  properties.  Thus 
Weibull  theory  is  based  mainly  on  extreme  value  statistics. 

Ideally,  statistical  fracture  theory  should  also  give  consideration  to 
fracture  mechanics  and  material  microstructure.  Progress  made  along  these 
lines  in  recent  years  provides  fresh  insight  into  Weibull' s theories.  The 
primary  objectives  of  the  present  paper  are  to  outline  a more  general  and 
fundamental  treatment  of  the  fracture  statistics  of  a solid  containing  randomly 
oriented  microcracks,  and  to  show  how  it  reduced  to  Weibull' s results  when 
certain  simplifying  assumptions  are  made. 
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FUNDA MENTAL  EQUATIONS  OF  WEAKEST  LINK  THEORY 

Let  it  be  assumed  that  c.  stressed  solid  can  fail  due  to  any  of  a number 
of  independent  and  mutually  exclusive  mechanisms  or  causes,  each  having 
an  infinitesimal  probability  of  failure  . The  probability  that  the  i'th 

mechanism  will  not  cause  failure  is  = 1 - ^ AP . The  overall 

probability  of  survival  is  the  product  of  the  individual  probabilities  of 
survival,  i.e. 


P 

s 


= n 

i 


2 U exp 
i 


? (APr). 


(1) 


The  sum  of  the  individual  probabilities  of  failure  appearing  in  the  final 
equality  above  was  called  oy  We’bull  the  "risk  of  rupture"  and  was  given 
the  symbol  B.  Note  that  the  "risk  of  rupture"  so  defined  can  be  arbitrarily 
large,  but  regardless  of  its  value  0 s Pg  Si,  as  of  course  it  must  on 
physical  grounds. 

Let  us  assume  that  the  potential  causes  of  failure  are  many  individual 
cracks.  For  purposes  of  analysis  it  is  convenient  to  group  the  cracks  accord- 
ing to  location,  the  applied  stress  state,  and  crack  critical  stress.  We 
assume  that  the  stress  state  varies  slowly  so  that  within  a volume  element  A V 
alL  cracks  will  be  subject  to  essentially  the  same  macroscopic  stress.  We  also 
assume  that  the  material  is  macroscopically  homogeneous,  so  that  a function 
N (cr)  depending  only  on  the  material  can  be  defined  as  the  number  of  cracks 
per  unit  volume  having  a critical  stress  equal  to  or  less  than  <7C.  When 
N(ac)<l,  we  can  regard  It  as  the  probability  that  a crack  of  critical  stress 
50c  will  be  found  in  a unit  volume.  The  critical  stress  of  a crack  is  defined 
as  the  remote  stress  which  will  cause  fracture  when  applied  normal  to  the  crack 
plane.  The  probability  that  a crack  having  a critical  stress  in  the  range 
oc  to  or  + dOcr  exists  in  volume  element  AV  is  then  Av[dN (ac) /dac]dac . 
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If  such  a crack  is  present,  the  probability  that  it  will  fracture  depends 
on  its  orientation,  the  stress  state,  and  the  fracture  criterion.  We  assume 
that  there  is  a solid  angle  ft  such  that  a crack  will  fracture  if  and  only  if  its 
normal  lies  within  ft.  This  means  that  if  the  normal  lies  within  ft,  C^  > O^, 
where  e is  the  effective  stress  corresponding  to  the  fracture  criterion 
selected.  If  the  cracks  are  randomly  oriented,  the  probability  that  a crack 
will  fracture  under  the  applied  stress  T is  ft  (S,  c )/ 4ir. 

The  probability  of  failure  due  to  a crack  in  the  critical  stress  range 
dc  located  in  volume  element  AV  is  the  product  of  the  above  probabilities, 
i.  e.  , 


= AV 


dN(o 
c 

dc 


n<C’°c> 

4tt 


(2) 


Substituting  equation  (?.)  into  equation  (I  ) and  changing  sums  into  integrals, 
we  obtain 


P 

s 


P, 


exp 


■fAVf 


dc 


dN 

dc 


4ir 


(3) 


Now',  if  w e assume  that  the  flaws  are  isotropic  as,  for  instance,  spheri- 
cally distributed  local  impurities  would  be,  equation  (3)  becomes  much  simpler.  In 
simple  tension  C,  ft/4tT  becomes  0 or  1 depending  on  whether  c is  less  or 

greater  than  c . Over  the  range  0 < c <C,  ft/4ir=l,  so 
f c c 


S 

/ 


. dN  ft  _ M . , 
dc  j— — -7—  = N.  (c 
c do  4tt  i 

c 


where  N.  is  the  density  of  isotropic  flaws.  Thus, 


= 1 - exp 


[-/ 


dV  n(C  ) 


(4) 


(5) 


where  we  have  changed  N.iC)  to  n(C  ) to  agree  with  Weibull's  notation. 
Equation  (5)  is  identical  to  that  employed  by  Weibull  in  Theory  A. 
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Accordingly,  the  simplest  interpretation  of  Theory  A is  to  say  it  is  based 
on  the  assumption  that  the  flaws  are  isotropic. 

There  is  another  possible  interpretation  of  equation  (5).  This  is  to 
permit  the  flaws  to  be  cracks  but  take  n(cr)  to  be  the  number  of  cracks  per 
unit  volume  that  will  be  fractured  by  a uniaxial  stress  a applied  to  the 
specimen.  Some  cracks  with  an  intrinsic  strength  considerably  less  than 
c will  survive  because  their  crack  planes  are  parallel  or  nearly  parallel 
to  the  tensile  axis.  Thus  n(CT)  is  smaller  than  our  N(ff),  which  is  defined 
to  be  the  number  of  cracks  per  unit  volume  having  an  intrinsic  strength 
equal  to  or  less  than  a.  We  can  also  regard  N(0)  as  the  number  of  cracks 
per  unit  volume  that  would  be  fractured  by  an  equitriaxial  tension  equal 
to  0. 

Returning  now  to  the  assumption  that  flaws  are  cracks,  the  probability 
of  survival  can  be  found  using  equation  (3)  provided  we  know  N(C  ) and  u. 
The  latter  is  a function  of  stress  state,  the  critical  stress  of  the  crack 
under  consideration,  and  the  assumptions  made  concerning  the  form  of  the 
effective  stress  a . In  simple  tension  and  equibiaxial  tension,  analytical 
expressions  can  be  found  for  In  the  general  case,  we  find  Cl  by  inte- 

grating dfiover  the  range  in  which  a ^ > c . One  way  of  accomplishing 
this  is  to  integrate  over  the  entire  angular  range  but  include  an  operator 
which  goes  to  zero  when  0 ^ < ff^.  Thus, 

P exp 

s 


4u 


Iff 


d V dfklff  IRo  , a ) , 
c e c da 


dN 


(6) 


c J 


where  H(o  ,<?c.)  = 1 when  > <if 

0 when  a < a 
e c 


Wt:  now  carry  out  the  integral  over  with  the  result 


P 


exp 


dVdCl  N(<7  ) 
e 


(7) 
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If  we  assume  the  effective  stress  is  simply  the  component  of 
applied  stress  normal  to  a crack  plane  0 , equation  (7)  is  essentially 
the  same  as  the  basic  equation  of  We  bull's  Theory  B,  and  equation  (3) 
becomes  equivalent  to  the  basic  equation  in  the  paper  by  Batdorf  and 

Q 

Crose.  This  means  that  Weibull  Theory  B is  valid  under  the  assump- 
tions of  Batdorf  and  Crose,  i.e.  under  the  approximations  that  cracks 
are  non-interacting,  that  they  are  stable  until  such  a stress  is  reached 
that  they  fail  catastrophically,  and  that  they  are  shear-insensitive, 
i.e.  the  effective  stress  is  the  normal  stress  C^. 

Other  assumptions  concerning  the  fracture  criterion  can,  of  course, 

9-13 

be  made.  A number  of  authors  have  assumed  that  fracture  occurs 

when  the  tensile  stress  at  some  point  on  the  surface  of  the  crack  cavity 
exceeds  tire  intrinsic  strength  of  the  material.  The  effective  stress  in 
this  case  depends  on  the  crack  shape  and  may  also  depend  on  Poisson's 
ratio.  This  approach  permits  evaluation  of  the  contribution  of  com- 
pressive stresses  to  fracture.  Other  authors  have  assumed  that  fracture 

. . , , 14-17 

occurs  when  the  strain  energy  release  rate  reaches  a critical  value. 

The  corresponding  effective  stress  is  readily  determined  for  in-plane 

crack  extension,  but  differences  of  opinion  exist  on  how  to  treat  the  more 

realistic  case  of  out-of-plane  crack  extension.  The  fracture  statistics 

in  polyaxial  tension  have  been  worked  out  for  a number  of  different 

~ 18 
fracture  criteria  for  shear- sensitive  cracks  by  Batdorf  and  Heinisch. 

Equations  (3)  and  (7)  are  alternative  and  completely  equivalent 
formulations  of  the  same  theory.  The  physical  justification  of  equation 
(3)  is  more  readily  apparent  than  that  of  equation  (7),  but  equation  (7) 
is  preferable  from  a computational  point  of  view. 
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FLAW  DENSITY  IN  THEORY  B 

We  have  already  discussed  flaw  density  n(CT)  in  fheory  A.  Weibull  did 
not  explicitly  introduce  flaw  density  into  Theory  B,  but  he  did  the  next 
thing  to  it.  In  Theory  A he  assumed 


Pj  = 1 - exp  r - J" n(cr)  dV 


n(cr)  = ka 

In  Theory  B,  he  assumed 


where 


Pf  = 1 - exp  j ~Jjni(an)  dVdft] 


, v , 111 

m(a)  = kia 

and  the  integral  extends  over  a half-sphere. 
Weibull  then  showed  that  for  simple  tension 

k i = (2m  + l)k/(2ir) 
or 

n i (a)  = (2m  + l)n(a)/(2rr) 


We  next  seek  a quantitative  relation  between  N(c)  and  n(a  ) which  we 
can  obtain  from  equation  (4)  when  we  know  fi.  Now  0 depends  on  the 
fracture  criterion  assumed.  If  we  assume  shear-insensitive  cracks, 


then  <3 


a , and  for  this  case  it  has  been  shown  that  in  uniaxial  tension 
n 


1 - V 


Inserting  equation  (14)  into  equation  (4)  and  integrating,  we  obtain 

N(a)  = (2m  + 1)  N.(o)  = (2m  + l)n(a)  . (15) 

Comparing  equations  (13)  and  (15),  we  see  that  Weibull's  n^  (o)  differs  only 

by  the  contant  factor  2tt  from  the  N(o)  in  the  present  paper.  He  did  not  give 

it  a physical  interpretation,  but  n^ (n)  can  be  regarded  as  the  number  of  cracks 

having  a critical  stress  less  than  a that  are  contained  in  1/(2tt)  units  of 

volume.  This  furnishes  additional  evidence  of  the  formal  equivalence  between 

8 

Weibull's  Theory  B and  the  theory  of  the  Batdorf  and  Crose  paper. 
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MATHEMATICAL  EXPRESSIONS  FOR  CRACK  DENSITY  FUNCTION 

More  often  than  not,  V\  eibull  assumed  that  n(ff)  takes  the  form  of  a 
simple  power  law  as  in  equation  (9).  This  leads,  in  the  case  of  a uniformly 
loaded  specimen,  to  the  result 


1 - 


exp 


V n (a)  = 1 - exp  - V ka" 


(16) 


In  practice,  a number  of  specimens  may  he  tested  to  obtain  Pf(cT),  and  k and 

m are  determined  by  a least  squares  fit  to  the  data.  The  probability  of 

failure  under  non-uniform  or  polyaxial  loading  is  then  found  by  inserting 

n(o)  = k a into  Eq. (8)  or  Eq.  (10).  The  result  of  integration  is  an 

equation  of  the  same  form  as  Eq . (16),  in  which  m is  unchanged  but  k is 
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modified  by  some  numberical  factor.  Batdorf  and  Heinisch  showed  that  this 
result  holds  not  only  for  the  shear-insensitive  flaws  considered  by  Weibull, 
but  also  for  shear-sensitive  cracks. 


The  use  of  a power  law  for  n(c ) implies  that  the  weakest  possible  flaw 
has  a strength  of  zero.  Weibull  also  considered  the  possibility  that  the 
strength  of  the  weakest  flaw  might  not  be  zero  but  . In  this  case  he 
suggested  that  failure  probability  in  simple  tension  should  be  expressed  by  a 
3 -parameter  relationship  in  the  form 


v rn  n 

r 

Pf  - 1 - CXP ] 

[-•(•.I  i 

= 1 - exp 

- V k(cr  -e  )m 

U 

This  introduces  no  serious  complications  for  Theory  A,  and  the  probabilities 
of  failure  under  variable  uniaxial  stresses  can  be  given  in  equations  of 
analogous  form.  However,  there  are  difficulties  in  the  case  of  Theory  B,  i.e. 
biaxial  and  triaxial  stress  states.  Basically  the  problem  arises  from  the 
fact  that,  as  Vardar  and  Finnie  have  pointed  out^,  the  integral  of  (c^  - O^)™ 

over  12  called  for  in  Eq.(7)  does  not  result  in  a constant  times  (o  - O )m. 

19  U 

Using  Theory  B,  Dukes  carried  out  extensive  calculations  for  various  values 

of  stress  ratios  and  values  of  a /a  and  expressed  the  results  in  graphical 
20  U ’ 

form.  Evans  has  carried  out  analogous  calculations  for  pennyshaped  cracks 
for  equibiaxial  and  equi triaxial  stress  states,  but  using  maximum  strain  energy 
release  rate  as  the  fracture  criterion. 
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If  the  cumulative  probability  of  failure  in  equation  (15)  is  differentiated 
to  get  the  probability  density,  the  result  within  the  range  0.01  < < 0.99  is 

a somewhat  skewed  Gauss  curve,  with  a skewness  determined  by  the  value 


of  m.  Consequently,  with  the  proper  choice  of  m,  the  equation  can  be 

used  to  describe  not  only  fracture  statistics  but  also  many  other  phenomena 

quite  accurately,  as  Weibull  pointed  out  in  1951.  It  is  not  clear  that 

Weibull  attributed  any  fundamental  physical  significance  to  the  form  of 

°2 

equation  (1  7),  but  there  are  some  who  have.  Epstein“  was  among  the 
first  to  point  out  that  this  equation  coincides  with  one  of  the  three 

23 

asymptotic  forms  of  extreme  value  statistics.  It  has  been  suggested 
that  all  statistical  fracture  theories  must  reduce  to  one  of  these  forms 
in  the  limit  of  large  specimens  (i.e.  V -*  ®).  We  shall  show  in  the 

next  section  that  this  is  not  the  case.  First,  however,  we  enquire  why  it  has 
appeared  plausible. 


Let  us  assume  that  ti(cr)  can  be  expanded  about  in  the  rather  general  form 


n (o)  = (o  - a )m  f.  a.  (a  - o )x  (18) 

u 1=0  1 u 

where  a = k % o,  and  m is  a positive  number, not  necessarily  an  integer.  We 
now  ask  how  P(.(a)  behaves  in  a region  of  interest  such  as  0.01  < <0.99  when 

v'  becomes  very  large.  From  Eq.(16)  it  follows  that  0.01  < Vn(c)  < 4.61.  Clearly 
as  V m,  n(o)  -*■  o,  so  that  (a  - a ) * o and  n (0)  k(o  - O^)  . 

Thus  from  Eq.(16a) 


->■  1 


exp 


V k (a  - a ) 
u 


(19) 


Now  equation  (19)  is  the  third  asymptotic  form  of  extreme  value 
statistics,  it  should  be  noted,  however,  that  this  equation  was  obtained 
as  a result  of  the  assumption  that  N (a^)  could  be  expanded  in  a gener- 
alized Taylor  series,  as  in  equation  (18).  While  equation  (18)  is 
plausible,  it  is  not  necessarily  a valid  description  of  crack  density  in 
real  materials.  In  fact,  as  we  shall  show  next,  it  probably  is  not. 
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McClintock  was  the  first  to  show  that  use  of  a physically  plausible 
model  for  crack  size  can  lead  to  a probability  of  failure  that  does  not 
reduce  to  any  of  the  asymptotic  forms  of  extreme  value  theory.  This  was 
accomplished  with  the  aid  of  a two-dimensional  material  model.  McClintock 
assumed  that  the  probability  that  any  two  adjacent  grains  are  not  bonded  is 
some  small  number  p,  and  that  a long  crack  is  simply  a random  aggregation 
of  unbonded  pairs  of  grains.  He  then  showed  that  the  probability  that  a 
crack  is  longer  than  x is  equal  to  exp  [-  x/X.]  . From  the  statistics 

of  crack  size  and  equations  of  fracture  mechanics,  he  derived  the  statistics 
of  critical  stress,  and  from  this  the  statistics  of  failure.  Earlier,  Fisher 
and  Holloman  had  worked  out  a statistical  theory  of  fracture  for  materials 
with  randomly  oriented  pennyshaped  cracks,  using  the  judgmentally  based 
assumption  that  the  probability  that  a crack  has  a radius  greater  than  r 
is  equal  to  exp  [-  r/\J.  The  resulting  theories  were  basically  very 

similar,  but  the  similarity  is  obscured  by  the  fact  that  Fisher  and  Holloman 
retained  the  binomial  product  form  in  their  calculations  rather  than  con- 
verting to  exponentials  [see  equation  (!)]•  Thus,  the  relation  of  their  work 

to  Weibull's,  and  the  implications  of  their  work  with  respect  to  asymptotic 

26 

forms,  escaped  general  attention.  Batdorf  applied  the  McClintock  hypo- 
thesis concerning  crack  origin  to  obtain  crack  size  statistics  for  penny- 
shaped cracks,  and  showed  that  it  leads  to  a probability  of  having  a crack  of 
radius  greater  than  r that  is  equal  to  exp  j^-  (r/a)^  j.  The  resulting 

statistics  for  N(c  ) take  the  form 
c 

N(Oc)  = NQ  expT-A/j^J  , (20) 

where  N^  is  the  total  number  of  cracks  per  unit  volume  and  A is  a somewhat 
complicated  function  of  p,  grain  diameter,  and  fracture  L ghness.  We  note 
that  Ntc^)  0 only  when  0,  but  that  N(c_)  cannot  be  expanded  as  a Maclaurin 
series  about  0 because  at  that  point  all  its  derivatives  vanish.  Thus,  we 

would  not  expect  P^.  (a ),  as  given  by  this  crack  density  function,  to  be  asymptoti- 
cally stable  and,  in  fac  t,  it  is  not.  The  theory  leads  to  the  following  result  for 


k- J 
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small  failure  probabilities 


- B j cr4  e 


(20) 


The  deviations  between  this  theory  and  Weibull  theory  are 
not  as  large  as  one  might  expect  from  the  dramatically  different  equations 
for  the  density  function  and  the  probability  of  failure.  Probability  of 
failure  curves  have  been  constructed  for  various  values  of  the  parameters 
(N  , p,  Kjc  , etc.  ),  and  it  has  been  shown  that  use  of  Weibull's  3-parameter 
equation  leads  to  a good  fit  over  a limited  range  of  failure  probability  (such 
as  0.  01  s Pf  < 0.  99).  However,  the  Weibull  parameters  and  m giving 
the  best  fit  are  weak  functions  of  specimen  volume  instead  of  being  volume- 
independent  as  in  Weibull  theory.  In  addition,  the  two  classes  of  theory 

differ  markedly  in  their  extrapolation  to  very  low  values  of  P^(a). 

It  is  not  suggested,  of  course,  that  the  theory  for  intergranular  cracks 
just  described  is  an  accurate  portrayal  of  physical  reality.  Several 
approximations  are  made,  such  as  the  assumptions  that  neighboring 
grains  are  either  perfectly  bonded  or  completed  unbonded,  that  cracks 
are  non-interacting  and  pennyshaped,  etc.  Future  work  will  undoubtedly 
lead  to  improvements  in  the  model.  The  point  to  be  emphasized,  however, 
is  that  an  improved  theory  can  be  expected  to  differ  from  Weibull  theory  in 
much  the  same  manner  as  the  theory  just  outlined  does.  It  seems  likely 
that  physical  reality  differs  from  Weibull  theory  in  a somewhat  similar 
way,  but  possibly  by  a small  enough  amount  to  make  it  difficult  to  establish 
these  differences  experimentally  with  complete  confidence. 
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CONCLUDING  REMARKS 

Weibull  was  the  first  to  apply  the  weakest  link  concept  to  brittle  frac- 
ture. He  showed  that  any  weakest  link  theory  must  result  in  a failure 
probability  for  a body  uniformly  stressed  in  simple  tension  in  the  form 

= 1 - exp  [-  V n(ff ) ] , (23) 

where  n(cr)  is  a positive  and  monotonic,  but  otherwise  arbitrary,  function 
of  the  applied  stress.  He  proposed  two  computationally  convenient  repre- 
sentations for  n(C  ).  one  involving  two  parameters  and  one  involving  three. 
These  are  so  widely  used  that  often  experimentalists  present  their  data 
simply  by  stating  the  values  of  the  Weibull  parameters  leading  to  the  best 
fit. 

Weibull  also  showed  how  to  apply  his  theory  to  predict  the  probability 
of  fracture  in  arbitrarily  stressed  bodies  when  the  failure  statistics  for 
simple  tension  are  known.  His  work  was  primarily  founded  on  the  consider- 
ations of  extreme  value  theory,  with  no  explicit  attention  to  either  fracture 
mechanics  or  material  microstructure. 

In  the  case  of  uniaxial  stress  states,  Weibull  assumed  that  only  tension 
contributes  to  fracture.  As  a result,  it  was  not  necessary  to  adopt  a 
fracture  criterion.  This  theory  is  basically  valid  for  tensile  type  fractures, 
subject  however  to  certain  limitations.  Use  of  either 

the  2-  or  the  3-parameter  Weibull  representation  of  fracture  statistics  in  a 
stress  range  for  which  test  data  are  not  available  is  unsafe.  This  means  we 
have  no  justification  for  extrapolation  to  very  low  probability  of  failure  or  to 
specimens  very  much  larger  or  smaller  than  those  tested.  The  hazard 
involved  in  such  extrapolations  is  made  especially  clear  by  the  results  of 
the  theory  for  the  fracture  statistics  of  intergranular  cracks  discussed  in  the 
previous  section. 

To  deduce  the  polyaxial  stress  fracture  statistics  from  uniaxial  data, 
it  is  necessary  to  specify  a fracture  criterion.  Weibull  did  not  do  this 
explicitly,  but  his  recommended  procedure  implies  the  assumption  that 
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fracture  occurs  when  the  component  of  tensile  stress  normal  to  a crack 
plane  exceeds  some  critical  value. 

Recent  improvements  in  weakest  link  theory  include  incorporation  of 
more  accurate  fracture  criteria  and  explorations  of  the  implications  of 
fracture  mechanics  and  material  microstructure.  While  these  develop- 
ments are  important  theoretical  advances,  it  is  too  early  to  assess  fully 
how  much  they  will  improve  fracture  statistics  from  a practical  point  of 
view.  In  any  event,  it  remains  an  impressive  tribute  to  Weibull's  technical 
grasp  and  physical  insight  that  the  world  has  required  several  decades  to 
discover  how  to  advance  fracture  statistics  significantly  beyond  the  point 
he  had  already  reached  in  1939. 
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